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Abstract 

By using the Furnstahl, Serot and Tang's model, the effect of density depen- 
dence of the effective nucleon-nucleon-p-meson (NNp) coupling on the liquid- 
gas phase transition in hot asymmetric nuclear matter is investigated. A limit 
pressure piim has been found. We found that the liquid-gas phase transition 
cannot take place \ip > p\\m- The binodal surface for density dependent NNp 
coupling situation is addressed. 
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It is generally recognized that the liquid-gas (LG) phase transition of one component 
system is of first order. The chemical potential continues at the phase transition point but 
its first order derivatives, namely, entropy and volume, are discontinuous. But for a multi- 
components or multi-conserved charges system, as was pointed out by Miiller and Serot [|l| , 
because of the greater dimensionality of the binodal surface, the LG phase transition can 
be of second order, i.e., the entropy continues but the second order derivatives of chemical 
potential (for example, capacity) are discontinuous. An asymmetric nuclear matter has two 
components of proton and neutron, and two conserved charges of baryon number and the 
third component of isospin, will undergo a continuous second order phase transition. 

Obviously, because of charge independence, the basic difference between proton and 
neutron be isospin. The isospin dependent interactions of nucleon-nucleon-isovector mesons 
play the key role to address the LG phase transition. As was pointed out by our previous 
papers IQ-Q, if one employed the isospin independent model, for example, Welacka model 
^ or Zimanyi-Moszkowski model [Q, to investigate asymmetric nuclear matter, a lot of 
difficulties, e.g. Coulomb instability and negative asymmetric parameter in the vapor phase, 
will emerge. To overcome these difficulties, the isospin vector p-meson must be introduced. 
It can be shown that the chemical potentials of proton and neutron may depend on the 
third component of isospin when NNp interaction exists. A model without isospin vector 
p-meson, or even if it has p-meson, but the chemical potentials of the proton or neutron are 
still independent of the NNp interaction because the third component I3 of isospin be zero 
such as in a symmetric nuclear matter, the LG phase transition is still of first order. 

In fact, the chemical potentials of proton and neutron not only depend on I3 but also on 
the effective NNp coupling gp. Then the effective coupling gp is also essential for studying the 
LG phase transition because the chemical potentials determine the binodal surface directly. 
In our previous papers we have shown that the effective couplings depend on the 

density and temperature in a hot and dense nuclear matter. By using the Thermo Field 
Dynamics to calculate the three-lines vertices Feynman diagrams of NNvr, NNa , NNo; 
and NNp interactions, we have found that the effective couplings of g^^, g^, g^ and gp all 
decrease as the nucleon density increases. In an asymmetric nuclear matter, one can easily 
prove that the chemical potentials of nucleons have a term which is proportional to g^ and I3. 
This term has opposite signs for proton and neutron due to their different third component 
of isospin. Obviously, if gp depends on density, this term will change and then the chemical 
potentials of proton and neutron, as well as the binodal surface of LG phase transition will 
also be changed. The objective of the present paper is to investigate the effect of the density 
dependence of gp on LG phase transition. We will prove that if gp is a decreasing function of 
density, a limit pressure will occur, when p > pum, the coexistenced equations have no 
solution and the LG phase transition cannot be existed. The chemical potential of neutron 
will become a monotoneous functioin of asymmetry a in this case. 

To illustrate our result, we employ a model suggested by Furnstahl, Serot and Tang (FST) 
|rT|-|l4| recently. This model is an extension of quantum hadrodynamics and has been proven 
to be successful to explain many experimental properties of both nuclear matter and the 
finite nuclei in mean field approximation. The Lagrangian density of FST model under mean 
field approximation is 
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where g^, gp are, respectively, the couphngs of hght scalar meson a, vector meson u and 
isovector meson p fields to the nucleon, (j)Q, Vq, bo are the expectation values 0o =< >, 
< >= 6^oVo, < >= S^obo- The scalar fiuctuation field is related to S* by 5* (x) = 
— 4>{x) and Hg is given by = 4ifg/ (^^S'q), d the scalar dimension. By using the 



standard technique of statistical mechanics, we get the thermodynamic potential as |T5 
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where j3 = l/fc^T and the quantity z/j (i = n,p) is related to the usual chemical potential /ij 
by the equations 
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where P3 = Pp — Pn and the third component of isospin I3 = {Np — N^) /2 = V P3/2. The 
third term of the right hand side of Eq(0) and Eq(|^) depends on and g^. They have 
opposite signs and play the essential role to determine the LG phase transition. 

Having obtained the thermodynamic potential, all other thermodynamic quantities, for 
example, pressure p=— fi/V^, can be calculated. The two-phase coexistence equations are 



/^f(r,pf)=/ir(T,pr) 
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p'{t,p^)=p^{t,p\) 
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where subscripts of one phase L and V stand for liquid and vapor, respectively. The stability 
conditions are given by [fff] 
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( ^\ <0oT ( ^\ > f8) 
where JF is the density of free energy, a = {pn — Pp) / p the asymmetric parameter, and 

P = Pn + Pp- 

The numerical calculatioins have been done by adopting the parameters set T 1 of FST 



model |1 1-13] 'The parameters of set Tl are 



^,2 = 99.3, gl = 15A.5, (?J = 70.2 (9) 
m, = 509MeV, So = 90.6MeV 
C = 0.0402, r] = -0.496, d = 2.70 

1 /2 

Our results for gp = (70.2) =constant are shown in Fig.l and Fig.2 by solid curves. The 
Gibbs conditions and @ for phase equilibrium demand equal pressures and chemical po- 
tentials for two phase with different concentrations. The collection of all such pairs ai {T,p) 
and ^2 {T,p) form the binodal suface. In Fig.l, the chemical isobar vs. a curves at fixed 
temperature T=10MeV and p=0.100MeV(/m)~^ are labeled by A and A' for neutron and 
proton respectively. The two desired solutions form the edges of a rectangle and can be 
found by means of the geometrical construction shown in Fig.l [Q. The critical curves with 
T=10MeV and pcru = 0.165MeV(/m)~^ are shown in Fig.l by B and B' where the chemical 
potential curve arrive at a inflection point and the rectangle is degenerate to a line vertical 
to the a axis. The behaviour of the nuclear matter under isothemal compression, or in other 
words, the section of binodal surface at finite temperature T=10MeV are shown in Fig.2. 
The physical behaviour of this processes has been discussed by ref. ll|. Assume that the 
system is initially prepared with a = 0.6 (gas), during the compression, the two-phase region 
is encountered at point A, and the liquid phase emerges at point B. The gas phase evolves 
from A to D, while the liquid phase evolves from B to C. The system leaves the region 
of instability at point C, while the original gas phase is about to disappear. The critical 
point (CP), the point of equal concentration (EC) and the maximal asymmetry (MA) are 
indicated in Fig.2. 

Now we are in a position to extend our discussion to the case of gp density dependence. 
In fact, the effective masses of nucleons, effective masses or screening masses of mesons, and 
the effective coupings of NN-mesons are all dependent on density and temperature. We can 
sum the tadpole diagrams and the exchange diagrams for nucleon, the vacumm polarization 
diagrams for mesons and the three-lines vertex diagrams for effective couplings to get their 
density and temperature dependence 0J^-^|T5|-0 . But in order to illustrate our result more 



transparently, instead of the exact calculation of three-lines vertex, we introduce an ansatz 



9p = gp 



I- Ap + Bp^ (10) 



where A, B are two adjust parameters. The reason for our choice are: at first, the three 
lines vertex calculations are model dependent, but we hope that our investigation can be 
more general; secondly adjust the values of A, B can made g'^ be decreased or increased with 
density, then we can study the LG phase transition for two different cases. We can imagine 
that the Eq(|TUp is an density expansion of effective coupling g'^ at low density regions. 

The results for density dependent g'^ ansatz Eq(p!OD are shown in Fig.3, Fig. 4 and Fig.5. 
We see from Fig.3 and Fig.4 that the chemical potential of neutron /i„ increases rapidly 
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with density. It passes through an inflection point and becomes monotoneous when pressure 
increases. But the shape of Hp vs. a curves change slowly. Then we will get a limit pressure 
Piim, when p > piiin, the rectangle cannot be found and the coexistenced equations have no 
solution. The last rectangle in the chemical isobar vs. a curves for A=l, B=0, T=10MeV, 
and piim=0.130MeV(/m)~^ is shown in Fig. 3 by dashed lines, where ai=0.62 and 0^3=0. 75 
correpond to the maximum and the minimum of fin respectively. The pair ai=0.62 and 
a;2=0.84 form the end of the binodal surface, as shown in Fig.5. The curve for A=l, B=0, 
T=10 MeV, but p=0.145MeV(/m)"^ {p > p^^) is shown in Fig.4. We see that /i„ becomes 
monotoneous at this pressure, and no rectangle can be found. The relation between limit 
pressure and parameters A and B for a fixed temperature T=10 MeV is shown in Table 
1. We find from Table 1 when A and B increase, the effective coupling g'^ decreases and 
the limit pressure decreases. If A changes its sign to become negative, g'p will increase with 
density, and in this case, instead of /i„ , fip becomes monotoneous. The limit pressure is still 
existed but decrease when A and B decrease. 

The section of binodal surface for A=l, B=0, T=10MeV is shown in Fig.5. We see from 
Fig.5 that the curve will cut off at limit temperature pum clearly. The total asymmetric 
parameter a is divided into four regions, namely, [0, ai], [ai, 03], [03, 0^2] and [02, Q^max]- The 
physical behaviour of isothermal compression in different regions are different. We find: 

(1) If the system is initially prepared with < a < ai, the precess of isothermal com- 
pression is similar to that of the case with constant gp. It begin at gas phase, suffers a second 
order LG phase transition and ends at liquid phase. 

(2) If the initial a is located at the region 02 < « < o^max, the system enters and leaves 
the two-phase region on the same branch, so the system remain in the same gase phase. 
As was pointed out by Miiller and Serot this retrograde condensation is unique to the 
binary system and dose not occur in one-component systems. 

(3) Suppose that the initial a is located at the region ai < a < 0:3. Since ai and 
correspond to the maximum and minimum of we find {diJ,n/da)rp^ < in this region 
and the stability condition Eq. (§) will be destroied. The system begins at gas phase, enters 
a two-phase region and becomes instable at the limit pressure. 

(4) If the initial a is located at the region < a < 02, the behaviou of the system is 
similar to that of the case (3), except it will be ended to a stable phase at the limit pressure 
because the stability condition {djin/da)j.^ > is satisfied. 

In summary, we have shown that the density dependence of effective NNp coupling is 
important to the LG phase transition. A limit pressure pum has been found for a fixed 
temperature and the LG phase transition cannot take place in asymmetric nuclear matter 
provided p > • Of course, for a fixed pressure, we can also get a limit temperature. This 
conclusion is similar to that of the Coulomb instability of nuclei. The basic difference 
is that instead of finite nuclei, our conclusion has be found to be suitable for asymmetric 
nuclear matter. Finally, we would like to emphasize that the isospin is very important for 
the LG phase transition of asymmetric nuclear matter. 

The work was support in part by the National Natural Science Foundation of China 
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TABLES 



TABLE L relation between the limit pressure and adjust parameter A, B 



A 


1 


2 


3 


5 


-2 


-2 


-5 


B 





1 


1 


2 


1 


-1 


-2 




0.13 


0.115 


0.105 


0.90 


0.125 


0.125 


0.115 
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FIGURES 



FIG. 1. The chemical isobar as a function of a at fixed temperature T=10MeV and the 
geometrical construction used to obtain the asymmetry parameters in the two coexisting phases. 
B, B' correspond to the critical curves. 

FIG. 2. The section of binodal surface at T=10MeV, where CP, EC, and MA stand for critical 
point, equal concentratioin and maximal asymmetry respectively. 

FIG. 3. The chemical isobar as a function of a. for A=l, B=0, T=10 MeV and p=0.145 
MeV(/m)~^ 

FIG. 4. The chemical isobar as a function of a for A=l, B=0, T=10 MeV and p=0.13 
MeV(/m)~^ 

FIG. 5. The section of binodal surface for A=l, B=0, T=10MeV. ai, a2 correspond to the 
maximum and minimum of /x„ respectively. 
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Fig5 
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